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Abstract: We study the cosmological evolution of the universe when quintessence is 
modeled within supergravity, super symmetry is broken in a hidden sector, and we also 
include observable matter in a third independent sector. We find that the presence of 
hidden sector supersymmetry breaking leads to modifications of the quintessence potential. 
We focus on the coupling of the SUGRA quintessence model to the MSSM and investigate 
two possibilities. First one can preserve the form of the SUGRA potential provided the 
hidden sector dynamics is tuned. The currently available limits on the violations of the 
equivalence principle imply a universal bound on the vacuum expectation value of the 
quintessence field now, k^I'^Q <^ 1. On the other hand, the hidden sector fields may be 
stabilised leading to a minimum of the quintessence potential where the quintessence field 
acquires a mass of the order of the gravitino mass, large enough to circumvent possible 
gravitational problems. However, the cosmological evolution of the quintessence field is 
affected by the presence of the minimum of the potential. The quintessence field settles 
down at the bottom of the potential very early in the history of the universe. Both at the 
background and the perturbation levels, the subsequent effect of the quintessence field is 
undistinguishable from a pure cosmological constant. 
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1. Introduction 

There is a host of observational evidence in favor of the existence of a non-zero vacuum 
energy density of the Universe driving the acceleration of the expansion of the Universe [1- 
3] . The simplest explanation for this new era in the history of the Universe is the presence 
of a cosmological constant of extraordinarily small value, some 120 orders of magnitude 
lower than the Planck scale. Such a small value is particularly difficult to accommodate 
when dynamical effects such as the Quantum Chromo-Dynamics (QCD) and electroweak 
phase transitions or even Grand Unified Theory (GUT) scale physics are taken into account. 
This has prompted the possibility of using extra dimensional models such as self-tuning 
scenarios [4] or brane induced gravity models [5]. Unfortunately these alternatives have 
drawbacks such as hidden fine-tunings [6]. Within four dimensional physics, there is an 
experimental way of discovering whether the vacuum energy is a true constant of nature or 
the result of more complicated dynamical effects. Indeed very active experimental programs 
are dedicated to the analysis of the so-called equation of state of the dark energy sector 
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(the ratio between the pressure and the energy density). If the equation of state differs 
from —1 (and is greater than —1, otherwise see for instance Ref. [7]), then a plausible 
candidate for dark energy is quintessence [8-10], i.e. the dynamics of a scalar field rolling 
down a runaway potential. Of course, quintessence only accounts for the small and non- 
vanishing vacuum energy, it has nothing to say about the large cancellation of the overall 
cosmological constant. 

One of the most stringent requirements imposed on quintessence models is the existence 
of attractors, i.e. long time stable solutions of the equations of motion [9]. Indeed the 
presence of an attractor implies an insensitivity to initial conditions of the quintessence 
field for the vacuum energy now. For a large class of potentials, attractors leading to 
vacuum energy dominance exist provided their large field behavior is of the inverse power 
law type. Such potentials are known under the name of Ratra-Peebles potentials [8]. In 
these cases, the quintessence field reaches an attractor, only to leave it when dominating 
the energy content of the universe. This happens when the field is of the order of the 
Planck scale. 

This has drastic consequences on quintessence model building. Indeed, it requires a 
natural framework within which Planck scale physics is taking into account. Supcrgravity 
is a promising field theoretical arena where both particle physics and Planck scale physics 
can be described [11]. Models of quintessence in supcrgravity have been constructed [12-15] 
leading to interesting phenomenological consequences such as low values of the equation of 
state. In particular, the simplest model of quintessence in supergravity, often dubbed the 
SUGRA model in the literature, leads to the following potential 

F,uint(Q)=e'^«^/^+<^^, (1.1) 

with K = 87r/mpj and M^+" = A^^^m°2°/^ and where, in this equation, Q is canonically 
normalized. The quantity a is a free positive index and A is a dimensionless coupling 
constant and, in order to avoid any fine-tuning, we will always consider that A ~ 1. rric 
is the cut-off scale of the effective theory used in order to derive the SUGRA potential. 
Typically nic can be thought as the GUT scale but we will see that the Planck scale is 
also possible (and, sometimes, necessary). Finally, ^ is a vacuum expectation value (vev) 
of some other field present in the quintessence sector, see below for more details. As a 
specific example, ^ can be realized as a Fayet-Iloupoulous term arising from the Green- 
Schwarz anomaly cancellation mechanism [12,13]. The main feature of the above potential 
is that supergravity corrections have been exponentiated and appear in the prefactor. 
Phenomenologically, this potential has the nice feature that the equation of state uj = 
PQ I PQ can be closer to —1 than with the Ratra-Peebles potential when the field approaches 
its present value K^'^^Qnow ~ 1- Moreover, a small value for M can be avoided. Indeed, since 
the vev of the quintessence field is now of the order of the Planck mass, requiring that the 
quintessence energy density be of the order of the critical energy density pen ~ lQ~^'^^rn^^ 
implies that 

^ 10-122/(4+a) ^ ^12) 
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and, therefore, M can be a large scale (by particle physics standard) for very reasonable 
values of the index a. For instance, it is above the TeV scale for a > 4. This mechanism is 
reminiscent of a "see-saw" mechanism where a very small scale (the cosmological constant 
scale) is explained in terms of a large one (the scale M) and a very large one (the Planck 
scale ?TT.pj). Moreover, the scale ^ can have acceptable values. Prom the expression of the 
scale M, one obtains 



and for a > 11 and a cut-off rric of the order of the GUT scale, is above the TeV 
scale [12]. However, considering rric < rUp^ can also be viewed as problematic since, as 
already mentioned, the vev of the quintessence field tends to be of the order of the Planck 



Facing this issue, a natural choice could be nic = nip^. Then the above equation indicates 
that ^ needs to be fine-tuned as the same level as the cosmological constant even if the 
model remains different since the equation of state is redshift dependent. Moreover, in 
this situation, the choice of the Kahler potential strongly influences the shape of the scalar 
potential [13] since the SUGRA corrections are of order one. Hence, it is no longer possible 
to see the Kahler potential of Eq. ( |2.2D as a Taylor expansion but it should rather be 
considered as a specific choice. Another route is to argue that, in these circumstances, the 
quintessence energy density can easily remain less than m^, indicating that the theory can 
still be meaningful if the cut-off is interpreted as the cut-off to the energy density scales 
and not to the vev's of the fields. In any case, it is clear that there are some fine-tuning 
problems at this level even if it is arguably less acute than in the cosmological constant 
case. 

However, it is clear that the quintessence sector cannot be considered as disconnected 
from the particle physics standard model (or its extensions) and should be embedded in 
a more general structure. In Ref. [16], we have investigated the coupling between the 
quintessence sector and a hidden sector where supersymmetry is broken. A general formal- 
ism to calculate the corresponding implications was presented. On very general grounds, it 
was shown that, as a consequence of this coupling, the shape of the quintessence potential 
is changed and that the particle masses become dependent on the quintessence field which 
implies the presence of a fifth force, a violation of the equivalence principle and possibly, 
depending on the complexity of the model, a variation of the gauge couplings and of the 
proton to electron mass ratio. 

The main goal of the present article is to apply this general formalism to a concrete 
case, namely the SUGRA one, where the quintessence potential is described by Eq. (pTl]). 
We find that the coupling between the quintessence, observable and hidden sectors has 
tremendous consequences on the dynamics of the quintessence field. The first one is that 
the shape of the quintessence potential is modified. When the hidden sector fields have 
been stabilized, the potential acquires a minimum located at a very small value of the field 
(in comparison with the Planck mass). As a result, the model becomes equivalent to a 
pure cosmological constant as the quintessence field settles at the minimum of a potential 




(1.3) 



mass. In this case, it is difficult to control the shape of the Kahler potential, see Eq. (p.2|). 
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before Big Bang Nucleosynthesis (BBN). The mass of the field is also changed and becomes 
equal to the gravitino mass. The above conclusion is also true at the perturbative level 
as no growing mode is present despite the smallness of the Jeans length. On top of this, 
the energy scales in the quintessence sector have to be fine-tuned at a highest level than 
the cosmological constant itself. This makes the whole scenario very unappealing. On 
the other hand, preserving the shape of the SUGRA potential requires a fine-tuning of 
the hidden sector dynamics. If this is done, one may wonder whether the model is still 
phenomenologically acceptable. We show that, in this case, the scenario tends to be ruled 
out by local tests of gravity rather than by cosmological considerations. 

Our results have been obtained using the SUGRA model [12] of quintessence. Eff'ec- 
tively the only crucial ingredient is the fact that the potential in the quintessence sector 
reduces to the Ratra-Peebles form for small values of the quintessence field. Within this 
framework, i.e. a quintessence sector with a Ratra-Peebles potential at small values of the 
quintessence field coupled to a hidden sector gravitationally, our conclusions apply and are 
generic even if a complete scan of the parameter space ("1,3/2, ''Tii/21 "^q) has not yet been 
performed. 

The paper is arranged as follows. In the following section, i.e. section |2[ we discuss 
the coupling of SUGRA quintessence to the hidden sector of supersymmetry breaking and 
the observable sector. In particular we find that the generic quintessence potential has a 
minimum with a mass for the quintessence field of the order of the gravitino mass. On 
the other hand one can fine-tune the hidden sector dynamics to preserve the form of the 
SUGRA potential. In section ^, we investigate the implications for gravity experiments 
of this fine-tuning of the SUGRA potential. In section ^ we discuss the cosmological 
evolution of the generic case where the quintessence potential develops a minimum and 
show that the field settles at the minimum of the potential before nucleosynthesis. We 
then examine the cosmological perturbations around the minimum of the potential and 
show that there are no growing modes for quintessence perturbations. We then conclude 
and mention possible way-outs in section ^. 

2. Quintessence and Supersymmetry Breaking 
2.1 The Framework 

As described in the Introduction, the usual approach consists in taking into account the 
features of the quintessence sector only. However, the non-observation of super symmetric 
partners to the standard model particles implies that SUSY must be broken at a scale of 
the order of a TeV. One must also take into account the existence of an observable sector, 
modeled as the Minimal Super Symmetric Model (MSSM) or mSUGRA model [11], where 
the standard model particles live. The calculations based on this more realistic description 
are performed in Ref. [16] where the model is separated in three sectors. There is a 
quintessence sector as already presented. It couples gravitationally to a visible sector where 
the standard model particles live. The supersymmetry breaking sector is also interacting 
gravitationally with the other two sectors. At the level of the Kahler K and super potentials 
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W, it is a simple sum of the contributions from each sector. 



K = i^quint + i^hid + Kobs , W = Wquint + Whid + VFobs • (2.1) 

In the present article, we use the general results obtained in Ref. [16] and apply them to a 
specific class of quintessence model. 

The quintessence sector is chosen such that it gives the SUGRA version of the Ratra- 
Peebles potential. It has been shown in Ref. [12] that this can be obtained using 

i^quint = QQ^ + xx^ + yytl^^ + ^ Uxl + y^y^) + • • • , (2.2) 



rric 



a=l 



Wquint = XX^Y + ^ XaXlY^ + ■■■ , (2.3) 



a=l 



where Xa and Y, Y^ are superfields satisfying {X) = ^, (Xa) = and (Y) = (Ya) = at 
the GUT scale where the model is defined (of course, in this context, the GUT scale is just 



an illustration), see however the discussion after Eq. (1.3). This implies that (Wquint) = 
and guarantees the positivity of the potential in the quintessence sector. The quantities 
A, Xa are dimensionless coupling constants and p is a free index. As already mentioned, 
the scale rric is the GUT scale below which the theory under consideration is valid. This 
means that the theory is valid only for vevs that are much less than rric or energy densities 
less than according to the interpretation given to the cut-off. The influence of the 
Kahlerian corrections are a priori important and have been studied in Ref. [13]. In the 
above expression and in the rest of this paper the dots stand for the higher order terms 
suppressed by mc- In the "dark energy" sector, we collectively denote the fields by da = 
X,Y,Xa,Ya,Q. Let us notice that we have assumed that the Kahler potential in the 
quintessence sector is regular at the origin. This excludes the no-scale case which deserves 
a special treatment, see Ref. [17]. It is interesting to consider regular Kahler potentials since 
they naturally lead to inverse power law scalar potentials while no scale Kahler potentials 
tend to give exponential potentials with very different properties, see Ref. [17]. 

For the hidden sector, we follow Ref. [16]. We denote the fields in the hidden sector by 
Zi and assume that Xhid is regular for small values of the hidden fields and can be Taylor 
expanded. Without specifying the superpotential for the moment, this leads to 

i^hid = Yl ^i^i + ■ • • ' ^hid = W^id{z^) , (2.4) 

i 

where, as before, the dots denote the higher order terms that are not considered in this 
article. 

Finally, the fields in the matter sector are written <j)a- This sector is supposed to contain 
all the (super) fields that are observable (including the dark matter). As a consequence, 
following again Ref. [16], we take this sector to be the MSSM or the mSUGRA model [11], 
that is to say 

-f^obs = ^ 4>a(t>i H , Wohs = Xabc4>a4>b4>c + ^ X] fJ-ab4>a4>b , (2-5) 

a abc ab 
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with a supersymmetric mass matrix ^ab and Yukawa couplings Xabc- In order to completely 
specify the observable sector, it is also necessary to choose the supergravity gauge coupling 
functions f^. A priori, all the /^'s are Zi and da-dependent. If, indeed, these functions 
are not constant, then this implies variations of the coupling constants. In particular, this 
leads to a variation of the fine structure constant, see Ref. [16]. 

Let us now discuss the breaking of supcrsymmctry in more details. In the hidden 
sector, the supersymmetry breaking fields Zi take a vev determined by 

d.,V = (2.6) 

where V is the total potential obtained from the previous model. The presence of the 
quintessence field affects the dynamics of the hidden sector and the vev's of the hidden 
sector fields become a priori Q dependent. They are parameterized in a model independent 
way by the coefficients ai{Q) and Ci{Q) according to 

(^.)^.^ ^ «.(Q) , «(l^,idLin~M3(Q), kV2/^\ ^ c{Q)M,{Q) , (2.7) 

\ / mill 

where Mg is the supersymmetry scale. Notice that, if the cut-off of the theory is much less 
than mpj, then we expect Oj <C 1 since the vev's of {zi) must be at most of the order of 
rric- Notice also that, a priori, nothing can be said about the function Ci{Q). On the other 
hand, if the cut-off is close to the Planck mass, the dynamics of the hidden sector cannot 
be analyzed in a model independent way and we must rely on particular models for the 
hidden sector to go any further. This is not surprising (and is even to be expected) since the 
situation is in fact exactly similar to what happens in the standard case of the mSUGRA 
model where the hidden sector is often described by the Polonyi model and the hidden 
field is stabilized at a vev of the order of the Planck mass. However, in the standard case, 
our ignorance of the hidden sector is parameterized in terms of two numbers, Oj and Cj, or 
equivalently ^3/2 and mi/2- In presence of dark energy, the situation is more complicated 
since our ignorance of the hidden sector is now described by two free functions ai{Q) and 
CiiQ). 

Another quantity of interest is of course the gravitino mass. It is defined by the 
following expression 

m3/2 = (AcW^e'^^/2) , (2.8) 

where K and W are the total Kahler and super potentials (i.e. taking into account the 
three sectors). In the present context, the gravitino mass may depend on the vev of the 
quintessence field. Therefore, it is natural to write 

m3/2 = e'^'^i-"*/2+S' l^'I'/^Mg = e''-^'i-"*/2m°/2 , (2.9) 

where is the mass that the gravitino would have without the presence of the quintessence 
field. This quantity explicitly appears in the expression of the scalar potential. Let us no- 
tice that Wquint does no appear in the previous formula because one has (Wquint) = for 
the particular case of the SUGRA model considered here. 
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2.2 The Dark Energy Sector 

In this subsection, we study how the shape of the SUGRA model is changed by supersym- 
metry breaking. As discussed in Ref. [16], see Eq. (2.18) of that article, the new potential 
in the dark sector is given by 



quint 



quint 



K 



-l\dl,di3 SXquint ^Xquint 3 



quint 



quint 



dd 



■/3 




quint 



quint 



(2.10) 



where V^uint is the potential that one would have obtained by considering the dark energy 
sector alone and the F-term in the hidden sector which takes the form 



^1/2 



(Mg +K(Wquint))ai+MgQ 



(2.11) 



In the case of the model described by Eqs. ( |2.2D , the above expressions simplify a lot, as 
(VFquint) = 0. Instead of the usual shape of the SUGRA potential, V{Q) = e'"(^'+^'+^" 
X Af^+^PQ"^^, see Eq. ( |1.1[ ) (in the minimal approach, the fields Xa are not present which 
explains the absence of the term X^a^o the above mentioned formula), we now have 



Q2p ' V 3/2 

where the function T(Q) encodes our ignorance of the hidden sector and is given by 



(2.12) 



r\Q) 




a=l 



+ e- 



\ 2 



^[a^iQ) + CiiQ)] 



— m. 



h/2 



(2.13) 



and M^+^P = e^* '""'^A^i^^mc^. In fact the quintessence field is not correctly normalized 
in this expression. To obtain a correctly normalized field, it is sufficient to replace Q with 
Q/^/2. The new shape of the potential is still not fixed in the above expression and is 
only known when the functions ai{Q) and Ci{Q) are specified, i.e. when the hidden sector 
is known explicitly. 

We can envisage two different situations. They are distinguished by the equation of 
state wq of the quintessence sector when Q takes its present value in the history of the 
universe. First of all, let us assume that the equation of state is wq 7^ —1. This can only be 
achieved when the potential is of runaway type with an effective mass for the quintessence 
field mq ~ Hq the present Hubble rate. This situation can only be achieved when the 
functions ai{Q) and Ci{Q) are not constant and such that, despite the new terms which 
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modify the shape of V^^, the runaway shape of the potential is preserved. This means 
that the susy breaking fields Zi are not stabilized or, more precisely, that the fields follow a 
trajectory in the field space {zi, Q). In this case, as we discuss in the following, one should 
precisely evaluate how serious the gravitational problems are. Obviously, this cannot be 
done in detail unless the functions ai{Q) and Ci{Q) are known exactly. In the following we 
will emphasize the case where aj <C 1 and q is tuned to obtain a runaway potential. 

Another case of particular interest is when a satisfactory model of the hidden sector 
has been found and the z^'s are correctly stabilized, i.e. ai{Q) and Ci{Q) are independent 
of Q. In addition, if we assume that the fields in the hidden sector are stabilized at 
vev's compatible with the cut-off rUc of the theory, i.e. we assume that {zi) <^ nic, and if 
rric <C rripj, then the coefficients are very small (since rric <C mp^) but this does not imply 
anything about the coefficients q. On the other hand, if nic is not small in comparison to 
mpj, then nothing can be said about Oj. 

Let us now focus on the simplest model of supersymmetry breaking where there is 
only one field z with a flat Kahler potential and a constant superpotential, PFhid = [18]. 
Then, one can even justify the previous choice. Indeed, first of all, this immediately implies 
that Mg = Km^ and c = 0. Then, the total potential (for simplicity, here, we assume 
Xa = 0) reads 



KZZ 



a=l 



'(QQt)P 



(2.14) 



where the dots now indicate the part containing the observable sector terms. These terms 
do not play a role at high energy and therefore can be ignored in the present context. At 
the minimum, we have dV/dz^ = 0, that is to say 



KZZ 



Km6(QQt)P 



0. (2.15) 



The constraint coming from the smallness of the vacuum energy implies (this is just another 
manifestation of the fact that quintessence has nothing to say about the cosmological 
constant problem) 

e+jiii--. (2.16) 

^ — ' K 

and, therefore, the quantity in the squared bracket is positive. As a result, the only solution 
is 

(^)min = 0, (2.17) 

i.e. a = 0. Hence, we find that this simple model satisfies = q = and, moreover, 
gives a SUSY breaking scale that does not depend on the quintessence field. In this case 
the function T{Q) becomes a constant and the potential in Q admits a minimum. At this 
minimum, it is clear that the mass of the quintessence field is of the order of the gravitino 
mass, mq ~ m^j2 = 0(l)TeV. In more complex settings, one can envisage a case where 
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one of the previous results is relaxed, for instance where c{Q) becomes Q-dependent while 
the other functions describing the hidden sector remain constant. This is what will be done 
in the following. 

We have thus reached one important conclusion, namely that when the quintessence 
sector corresponds to the SUGRA potential and when the hidden sector is correctly stabi- 
lized then the quintessence field acquires a mass of the order of the gravitino mass. If the 
hidden sector is more complicated, the runaway shape can be preserved and it is interest- 
ing to see whether local gravity tests can constrain this type of models. The cosmological 
implications of the new potential ( ^.12 ) are worked out in detail in Sec. 



2.3 The Observable Sector 

After having described how the dark energy sector looks like for the SUGRA model, we 
now study the observable sector. On very general grounds, it was shown in Ref. [16] that 
the interaction between the dark energy sector and the observable sector implies a Yukawa 
like interaction of the form 

2 f Q 



m 



(2.18) 



where ^' is a fermionic field. Therefore, the above expression implies that the fermion 
masses become quintessence field dependent quantities. As we discuss in the following, 
this implies a series of interesting effects as the presence of a fifth force and/or a violation 
of the weak equivalence principle. Our goal in this subsection is to use the general results 
obtained in Ref. [16] and to apply them to the SUGRA model in order to compute explicitly 
the functions m{Q/m^^). 

Computing the masses of the fermions first requires to compute the soft terms. This 
was done in Ref. [16], see Eqs. (2.21)-(2.23) of that article, and the general result reads 



A 



abc 



B. 



ab 



m 



ab 



ddr- 



dd, 



quint f^^'^^quint 



ddi 



ddf; 



(2.19) 




quint 



ddi 



dd 



■/3 



-,kK 

quint I 



+ Mg ( kW, 



quint + /^W^quint 



s. 



ab ' 



(2.20) 



(2.21) 



'quint I' I' quint 

This is the general form of the soft terms, calculated at the GUT scale. Then, we have to 
specialize the above formulas to the dark sector described by Eqs. ( |2.2[ ). Again the fact 
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that (Wquint) = considerably simplifies the calculations. Straightforward manipulations 
lead to 



Aabc — -^abc^3/2^ 



kK. 



+ ^ ^ aiCi + ^ ( kQ"^ + Ki'^ 



a=l ^ 



(2.22) 



IJ'ablTT'S/2 
n 



-Squint g. 



ab 1 



(2.23) 
(2.24) 



where, again for simplicity, we have assumed = 0. In the following, one will neglect ^ in 
comparison with the vev of Q (as ^ has to be extremely small in order for the quintessence 
energy density to be of the order of the critical energy density today, see above) and one will 
also consider that ^1^=1 = This last relation appears in the model of the hidden 
sector discussed before, see Eq. ( 2.16|) , but is also natural in the present context when one 
treats the case where the potential keeps its runaway shape despite the appearance of the 
new terms coming from the hidden sector. This means that V^^ should vanish at infinity 
and, therefore, that the constant terms cancel in the term T^(Q). Again, this issue is 
linked to the cosmological constant problem. 

In order to obtain the masses of the fermions, one should use the following procedure. 
In the MSSM, they are two Higgs doublets and and the fermions either couple to 
the Higgs doublet "u" or "d". These couplings are different and, therefore, through the 
Higgs mechanism, the masses of the fermions either depend on the vev of the Higgs "u" or 
of the Higgs "d" . Exphcitly, one has [16] 



;^Fgi^quint/2+EJa«lV2 



md 



yjl + tan^ [5 
v{Q) + 

;^Fgi^quint/2 + E, I«.IV2 V 

x/l + tan2 f3 



;^Fgi^quint/2+EJa«lV2 



v{Q) 

tan (3 



+ 



tan2/3 



tan^/? 



(2.25) 



(2.26) 



where ^ are the Yukawa coupling which are, in the minimal setting considered here, 

independent of the quintessence field. In the previous equations, v is defined by \J v1^ + 
where Uu = f sin (3 and Vd = v cos (3 are the vevs of the Higgs and respectively. We 
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have expanded the two previous expressions in terms of 1/ tan/3. The exphcit expression 
of tan/3 reads [16] 



tan/3(Q) 



2tiB{Q) 



li-^/l- V52(g) 2|/x|2eS.kP+m2^JQ)+m|^^(Q) ).(2.27) 



Notice that there are two possibihties according to the sign in the above expression. In 
the following, we work in the limit of large tan /? and, therefore, take the largest value. 
In Eq. (p.27| ), m'jj (Q) and m'jj (Q) are the two loops renormalized Higgs masses given 
by [16,19] 



1 



^hAQ) =<,(Q)-0.36 1 + 



d'"" " \ ' tan2/3 
+ ( 0.28 - ) A{Q) + 2m^/2 



0.72 



tan2/3 



iQ) 



m. 



'3/2 



1 



0.15 

in 



+ 2 i"'i/2 




We see from the two above expressions that the quintessence field dependence of the 
fermions masses is controlled by two functions, A{Q) and B{Q). The explicit expres- 
sions of these soft terms, see Eqs. (2.22)-(2.24), allow us to extract the functions A[Q) and 



B{Q) by Aabc = e*^ 'i"'"M(Q)Aafec and Bab = e*^ 'i"'"'^-B((5)ea6. If we restrict our consider- 
ations to aj = 0, which is a case of interest since this is at the same time compatible with 
{zi) <C TTip and with a possible runaway potential (using a non trivial dependence of the 
coefficient q, see the discussions before), one arrives at 



A{Q) =mAi + 



B{Q) = mAi + 



(2.30) 



Notice that A and B follow the same universal relationship as in the mSUGRA model 
despite the presence of the quintessence field. Then, one deduces that tan/3(Q) in the 
SUGRA model can be expressed as 



tan/3(Q) ~ 



5^jf52KQ2^f53_K2g4 

(54 -I- b^KQ'^ 



1 + Wi + 



(<54 + <55kQ2)2 



(5i + h^Q^ + (53k2Q4)2 



(2.31) 



where the coefficients 82-, i^s, 64^ and ^5 can easily be evaluated in terms of the physical 
parameters characterizing the model from the previous equations, i.e. /i, m^j^ and ?t^5/2 
given at the GUT scale. As already mentioned, this result is valid both in the case where 
the hidden sector fields are stabilized with Oj = Cj = and when the hidden sector dynamics 
is tuned to reach Oj = and Q 7^ possibly leading to a runaway potential. In the latter 
the role of the gravitational tests is crucial in discriminating models. 
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The expression of the scale v{Q) can also be obtained from the minimization of the 
Higgs potential along the lines described in Ref. [16]. One obtains 



v{Q) = _=^e'^^--/%/ + m2 , + O ( ) , (2.32) 



|2 ^^2 



tan/3 



where as before, we have used Oj = 0. The value of v today is known and is v ~ 174 GeV. 
Therefore, recalling that m^o = {g^ + 5'^)^^^/2 with m^o ~ 91.6 GeV and expressing 
Eq. ( 2.32| ) at vanishing redshift allows us to determine the n parameter. Explicitly, one 
has 



\H\ = ^ im|oe-'^('?— ) - . (2.33) 

In the following, as already explained, we neglect ^ and use Y12^=i^a = 3//t. The above 
formula also depends on the value of the vev of the quintessence field today, Qnow This one 
should be determined by the cosmological evolution of the field in the potential V^j, and is 
therefore fixed once the parameters (i.e. ""^3^2' "^1/2 • • ■ ) ^^'^^ been chosen. However, 
in the runaway case, it is clear that one always have K^/^Qnow ~ 1 and, for definiteness, 
we will take k^/^Q 

now — 1 in the following, the final results being (almost) independent of 
the precise value of i^^^'^ Qnow 

Let us also remark that the above formula giving tan /? and v are only approximated 
formulas and, as announced above, valid only when terms like 1/tan^/? are negligible in 
Eqs. (g) and ( ^ ). Otherwise one would have to deal with a transcendental equa- 
tion. If necessary, this equation can always be solved numerically, but, in this article, we 
always consider the approximation where the various quantities of interest are expanded 
in l/tan/3. The corresponding evolution of tan/3 as given by Eq. ( |2.31| ) is represented 
in Fig. |l[ We see that tan/? does not change sign, which would have been problematic 
once the sign of ^ is fixed. In Fig. [l|, one check that the electroweak symmetry breaking 
conditions are indeed satisfied in the whole range for which tan /? is plotted. We mentioned 



before that this plot has been obtained by neglecting the terms 1/tan /? in Eqs. ( 2.2 



and (2.29). One can verify that, for our choice of parameters, this is a good approximation 



since tan/? > 5. In Fig. g, we have also represented the scale v{Q) given by Eq. (2.32) for 
the same values of the parameters as before. 



If the expressions of tan/3 and v{Q) are used in Eqs. ( 2.25 ) and (|2.26|) , this gives 



^u,d{Q)- To our knowledge, this is the first time that the Q dependence of the fermions 
masses is calculated in a precise model from first principles. The vev Vu{Q) is equal to 
v{Q) at leading order and the vev v^iQ) is represented in Fig. ^ 

Finally, let us remark that, very often in the literature, the function mu,d(Q) is just 
postulated see, for instance, Ref. [20] where rrin^diQ) exp (/?k^/^Q/2) , /? being a free 
parameter. We see that, even in our oversimplified model, the dependence can be much 
more complicated. 
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Figure 1: Evolution of tan/3 versus tlie vev of the quintessence field according to Eq. ( 2.27| ) or 



Eq. (2.31) 



3. Implications for Gravity Experiments 



3.1 Fifth Force Constraints 

In this section, we study the consequences of the fact that the fermions masses are now 
Q-dependent quantities. It is known that this can cause serious problems coming from 
gravitational experiments since this implies the presence of a fifth force and a violation of 
the weak equivalence principle, see the next subsection. In fact, this crucially depends on 
the mass of the quintessence field. If the mass of the quintessence field is larger than 10~^eV, 
then the gravitational constraints are always satisfied as the range of the force mediated 
by Q is less than one millimeter. We see that this occurs when the functions ai{Q) and 
Cj(Q) vanish or are constant. Indeed, we have shown that, in this case, the potential has 
a minimum and acquires a mass ~ ^-3/2 10~^eV. Therefore, we reach the conclusion 
that the SUGRA model, with a hidden sector such that the fields are correctly stabilized, 
is free from gravitational problems. As discussed in the next section, the problems rather 
originate from cosmological considerations. 

On the other hand, if the mass of the quintessence field is less than 10~^eV, the range of 
the quintessence field is large and generically there will be violations of the weak equivalence 
principle and a large fifth force. In the case of the SUGRA model, this requires non trivial 
functions ai{Q) and Ci{Q). In this situation, in order to avoid fifth force experiments such 
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as the recent Cassini spacecraft experiment, one must require that the Eddington (post- 
Newtonian) parameter I7 — 1| < 5 x 10~^, see Ref. [21]. The Eddington parameter 7 is 
related to the parameters Ou^d by 7 = 1 + a^^, where a^^^ can be expressed as 



"u,d(Q) = 



1 dlnmf^^(Q) 
dQ 



din 



dQ 



(3.1) 



Therefore, the coefficients a^.d can be determined explicitly from the formulas giving mu,d- 
The model is free from difficulties if ^ < 10~^. Clearly, the coefficients can be calculated 
explicitly as soon as the functions ai{Q) and Ci{Q) are known. Here, we perform this 
calculation for the choice = leaving the function Cj 7^ free to lead to a runaway 
potential as already discussed at length before. In addition, in order to deal with the 
simplest model, we consider that the scale Mg does not dependent on Q as indicated by 
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the model of SUSY breaking presented before. One obtains 



_ k-V^ dtan/3 K-^&u_ 

«u - -^dgK^uint + tan/3(l+tan2/3) dQ + ^dQ 

dgK^uint + ^ + O , (3.2) 



2 "^"4-- ' t; dQ ' ytan'P 

- '^^^^ fi K K"^/^tan/5dtan/3 dt; 

«d - ^c^Qi^quint - i + tan2/3 dQ + ^^dQ 

^-1/2 dtan/? ^-1/2 dt; / 1 
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where the derivative of the function tan (3{Q) can be expressed as 



dtan/3 
dQ 



+ '^"'^^V2luP + m2 +mM-Han/? ^ ^^^^^ 



\^ dQ dQ 
±2^(2|/x|2 + m|^ + m2jJ-' 



tan P 



l-Af,^B\Q)[2\i2\-' + mi^+mi^ 



dB{Q) 
dQ 



+ B{Q) 



dm^^ 



+ 



dmy^ 



dQ dQ 



dm?r dm?, 

H-u _|_ Sid 



dQ 



dQ 



tan /? 



dlnB(Q) 
dQ 



tan (5 . 



-1/2 

(3.4) 
(3.5) 



The last expression is valid at leading order and we use the explicit form of the functions 
A{Q) and B{Q) for the SUGRA model in order to evaluate the derivatives of the soft terms 



d^(Q) 



2m°y2 



kQ . 



dB{Q) 



m 







(3.6) 



dQ 3 dQ 3/2' 

since Mg is constant. Consequently, the derivatives of the Higgs masses can be expressed 
as 



dm^ 
dQ" 



-0.72 X 0.28 



d^(Q) 
dQ 



A{Q) + 2m° 



1/2 



dm 



dQ 



0. 



(3.7) 



the symbol "approximate" in the last two equations meaning that we have used the fact 
that the terms in 1/ tan^ /3 have been neglected in the expression of the above formulas. 
In the previous calculation, we have also used the fact that m^^g is constant. This means 
that we have assumed specific forms for the gauge functions , namely we have considered 
that they do not depend on Q and Zi but only on the dark sector fields Xa- Then, if we 
parameterize and the derivative of / = /„ as 



3 

K 



-1/2 



dl_ 



where the coefficients e„ and hn are of order one, one finds that 



m 





1/2 



(3.8) 



(3.9) 



with no dependence on Q and a model dependent prefactor of ^-2^2- 

As an example we have plotted in Fig. Q the gravitational coupling constants a^^d for 
a realistic situation where o,; = and the parameters m^/o = 150 GeV, = 580 GeV 



and the parameters rriy^ 
and, hence, = 986 GeV, that is to say m? /n and m? roughly speaking of the same order 



''3/2 ^^^^ '"1/2 

of magnitude as indicated by the previous calculation. We see that the limit Ud ~ 10~^'^ 
is reached for relatively large value of the quintessence field vev, or the order of k^/^ (Q) ~ 
10~^ — 10~^. This implies that the SUGRA model with a, = and c, 7^ to obtain a 
runaway potential with Qnow 

~ mpj is excluded unless the dependence of the masses on Q 



in (3T) involves a Q dependent Yukawa coupling compensating exactly the Q dependence 
of the Higgs vevs. Although this is not excluded, this is a functional fine-tuning which is 
hard to explain. 



-16- 



■d 



0.010 



0.008 



0.006 



0.004 



0.002 



ai=0 

m3/2°=150 GeV 
mi/2°=580 GeV 
/^=986 GeV 



-7 — B -4 —3 -? -1 

10 10 ^ 10 ^ 10 10 ^ 10 ^ 10 



10 







Figure 4: Evolution of the coefficients au (solid line) and aa (dashed line) versus the vev of the 
quintessence field. 



3.2 Violation of the Weak Equivalence Principle 

As explained in detail in Ref. [16], the fact that, in the MSSM, the fermions couple dif- 
ferently to the two Higgs doublets and if^ causes, in the presence of dark energy, a 
violation of the weak equivalence principle. This violation is quantified in terms of the ry^g 
parameter defined by [22-24] 



Aa 



2^ 



a, + Op 



(3.10) 



for two test bodies A and B in the gravitational background of a third one E. Current 
limits [25] give 77^3 = (+0.1 it 2.7 it 1.7) x 10"-^^. The r/^g parameter was computed in 
Ref. [16] where a general formula was derived. Applying this general result to the case of 
the SUGRA potential leads to 



d_ 
dQ 



a 



A 



qCD 



+ 



d 



dQ VA 



5' 



'■QCD 



(3.11) 



where N^^ (respectively N^) represents the number of neutrons in the atom A, the body 
A, by definition, being made of this type of atoms (respectively in the atom B) and 
(respectively Z^) represents the number of protons in the atom A (respectively in the 
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atom B). The coefficients is defined by = K'^/^dlnm^/dQ but for two pairs of 
test bodies, the ratio rj ab/Vbc independent of the background object E. The quantity 
^QCD ~ ISOMeV is the QCD scale. Finally the variation of the coefficients a' and 5' can 
be expressed as 

d f a' \ 1 1 K^/^ 

(^u + &d) (au"iu + ad^T-d) + TT-r aame , (3.12) 



dQ X^QCU J ^ ^QCD ^ -'^QCD 

d f 5' \ 1 I ^1/2 

{K - fed) (OuWlu - ad"ld) - «d"ie , (3.13) 



dQ VAgj^py 2 Aqpq 2 Aqi-,j3 

where feu + ^d ~ 6, feu-^d ~ 0.5, ct7Aq(.q ~ 3.8 x 10"^, 5' / K^^^ ~ 4.2 x IQ-"^, mu ~ 5MeV, 
rrid ~ lOMeV and me ~ 0.5MeV. Let us notice that the coefficient appears in front of 
the mass of the electron mg because, in the MSSM, the electron behaves as a "d" particle. 
If we compare the above equations with the ones of Ref. [16], one sees that there is no 
variations of the fine structure constant and/or of the gauge function /qcjj. This is simply 
because, in this article, we assume that they are constant. 

Let us recall that when the mass of the quintessence field is of the order of the gravitino 
mass, all the gravity tests, among which is the equivalence principle violation are trivially 
satisfied. However, it is interesting to compute the r^^g in the SUGRA model with the non 
trivial hidden sector. One obtains 

1 . ^ fN^ N^+Z^\ I , 

Vab ~ (0.084a, + 0.168ad) [ ^ " ^ ) + (-0.0069a, + 0.013ad) 

when Qi = and q 7^ to lead to a runaway potential. It is interesting to notice that 
the calculations of the equivalence principle violation is directly related to the supergravity 
Lagrangian and, therefore, originates from first principle. A good order of magnitude 
estimate of the parameter r/^^ is simply, see Fig. § 

Vab ~ "u > (3-15) 



where we have assumed a^, ~ a, ^ aa- In order to comply with the currently available 
experimental limits this implies 

K^/^Qnow < 10"^ . (3.16) 

The above number is obtained for rn^^,^ = 150 GeV and = 580 GeV, which implies 

/X = 985 GeV, i.e. the values used in Fig. ^. This means that constraints on the weak 
equivalence principle are able to rule out this version of the SUGRA model. In other 
words, if the hidden sector is chosen such that it leads to an interesting cosmological 
model, then local tests of gravity are able to kill the scenario (at least for these values of 
the parameters). 



-18- 



3.3 The Proton-Electron Mass Ratio 



Another consequence of coupling the MSSM to dark energy is that this will cause the 
proton to electron mass ratio, r = rrip/me, to vary with time, that is to say with the 
redshift. This is an important consequence since, experimentally, there is an indication for 
a possible variation, Ar/r ~ (2.0 it 0.6) x 10^^ [27]. The proton mass can be written as 



rrir 



Cqcd^qcd + burriu + b^nid + Cpa^^^ 



5.2 and CpOQ^^ ~ 



''^ -Squint/ 2 



(3.17) 
0.63MeV. 

(3.18) 



where mu and are the mass of the u and d quarks, C^^-^ 
Recalling that the electron behaves as a "d" particle, this leads to 

r = tan/3 + b^— + VI + tan^ l3— ^ — tttt — « 

where Au, Aj and Ag are respectively the Yukawa couplings of the quarks u, d and of the 
electron. We also assume, for simplicity, that the gauge functions are constant, see above. 
Using Eqs. (p35| ) and ( ^) , one can express the Yukawa couplings in terms of the fermion 
masses at present time, generically denoted "m"". This leads to 



mr, 



tan^(Q) 



rn-O tan/?(Qnow) 
tan/?(Q) 



+ b 



"^d , ^QCD^QCD + ^pO^QED 
'd — 



tan/?(Qnow) ^ ' 

As before, in order to be consistent, we have to expand the above equation in terms of 
1/ tan (3. Then, one obtains 



+ 



b ^ + 
rrv^ 



^QCD-^QCD ^P'^QED 

m2 



tan/3(Q) 
tan/3((5now) 



1 



(3.20) 



One can also Taylor expand the tangent function and use the expression of the derivative 
of the tangent at leading order. This gives 



Ar 

r 



1 



dm^^ 



dm? 



rl/2 



+ 



dQ dQ 



1 dlnS 



,1/2 



Ag. (3.21) 



For instance, if one uses our fiducial model with m^i^ 

has ^l = 985.7 GeV, m^j^ 971650.6 GeV^ and m^^ ~ 

for order of magnitude estimate, we have taken i^^^"^ Qnov^ 
runaway behavior. This leads to 

Ar 



150GeV and m'^^^ 
190431.3 GeV^ where we recall. 



580GeV, one 
3re we recall, 
1 which is consistent with a 



0.69k^/^AQ, 

This means that, in order to comply with the recent bound, one should have 

K^/^AQ < 10^^ , 



(3.22) 



(3.23) 



in the range z € [0,3]. One sees that this constraint is less strong than the one obtained 
from the equivalence principle, at least for this choice of parameters, but is also sufficient 
to rule out the model. 
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3.4 Cold Dark and Baryonic Energy Densities 

The coupling between the observable sector (radiation and matter) and the quintessence 
sector also affects the way the energy densities scale with the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) scale factor a. In a minimal setting, it turns out the quintessence field 
couples to matter and not to radiation since the gauge functions are chosen to be 
constant. This implies that the radiation energy density still behaves like 

/Orad oc ^ . (3.24) 

In particular, in this case, there is no variation of the fine-structure constant. On the 
contrary since the masses of the particles depend on kQ^, the matter density behaves as 

/Cmat = ^ ^aml d [l^Q^] , (3.25) 

a 

where (^(kQ^) is the mass of non relativistic species the expression of which has been ob- 
tained before. The quantity is the number of non-relativistic particles which is conserved 
ha + SHua = 0, where H = a/a is the Hubble parameter. This leads to 



a 



where is just a constant. Again, in the most general case, the scaling of the cold 
dark and baryonic energy densities depends on the functions ai{Q) and Ci{Q) and is com- 
pletely explicit in the case Oj = 0. Let us also notice that cosmological consequences of 
Eq. ( 3.26| ) have been studied in Ref. [20]. As one can see, the coupling between matter and 



quintessence induces a modification of the quintessence potential 

KffCQ) =V;,e(Q)+^cdm(Q)%^, (3.27) 

where A^^y^{Q) = rn^J^y^{Q) /m^-^y^{fi), ""^cdm being the mass of the dark matter particle, 
typically the lightest supersymmetric particle. When the effective potential admits a time- 
dependent minimum, the model is known as a chameleon model [28] and the cosmological 
evolution is changed. However, in the present context, this correction is negligible in 
comparison with the correction coming from the susy breaking, i.e. m'^i^Q'^, see Eq. ( 2.12| ). 



In fact, one must compare the derivatives of the two corrections since this is the derivative 
of the potential which appears in the Klein-Gordon equation. For the new term coming 
from the dark matter energy density, one has 

and for the susy breaking term, this is simply m^i^Q. Since rn^i^ ^ Pmat/'"^pi, one can 
ignore the correction coming from the cold dark matter energy density as announced pre- 
viously. This is essentially due to the susy breaking term as the field happens to settle 
down at a minimum where its vev is very small. The situation would drastically change if 
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the potential kept its runaway shape as could be the case, for instance, if the hidden sector 
were such that ai{Q) and Ci{Q) were non trivial. In the following, we do not consider this 
situation and, therefore, we neglect the correction coming from the non trivial dependence 
of the dark matter energy density. 



4. Implications for Cosmology 

4.1 Fixing the Free Parameters of tlie Potential 

In this section, we study the cosmological implications for the case = Ci = as we 
have just seen that this case cannot be ruled out by the gravity constraints. The previous 
quintessence potential with = Cj = is characterized by three free parameters, M, 
777-3^2 and T. As a consequence, the potential can be rewritten in terms of dimensionless 
quantities as 

T^quint(Q) = m^e^'/^ + AQ^ - B) , (4.1) 

where Q = n^^'^Q is dimensionless and (obviously, the A and B below have nothing to do 
with the soft terms calculated in the previous subsection) 

M)'- . (s.r«. ,4(s.,w.(!4.y(My'-', ,,,, 

B.(s.,-«/.(JL)Vil)-<^^"'. ,,3, 

The above potential has a minimum, Qmin- Let us assume that we are in a situation where 
Qmin ^ 1- In this case, an explicit expression of Qmin can be found. Since the value of the 
field at the minimum is small, the exponential factor will be close to one. If so, it is easy 
to show that 

Qmin (I) . (4.4) 

We see that the location of the minimum is now controlled by the gravitino mass through 
the constant A. It is interesting to compare this value with the minimum of the usual 
SUGRA potential which reads Qmin = Va and is controlled by the Planck mass only. 
As we show in the following, the value of the minimum in the case where supersymmetry 
breaking is taken into account is much smaller than the minimum of the SUGRA potential. 
If we define the reduced field ^ by g = Q/Qmin, then the potential becomes 

Vquint(g) = m^e'Siin'/VS (^q-a + _ _ (4 5) 

The new scale rh is of course directly related to the scale m. The relation reads = 
Qmhi'^^- '^^^ potential is still characterized by three parameters which are now rh, B and 
Qmin- The constant rh is chosen such that the quintessence energy density be approximately 
70% of the critical energy density today. If we assume that the field is, today, at its 
minimum, i.e. (? = 1, (this will be shown in the following) this leads to 

^ \ — ^1^— s'^- ^ ^(1)'^- ■ (^■^) 

1 + a/2 — B 
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Then, one can link the remaining parameters with the particle physics parameters. For the 
gravitino mass, one obtains 

(<A\^Q-2 Pen , . 

As an example, one can take rn^^^ ~ lOOGeV which gives Qmin ^ 2.6 X 10 ^ wll6r6 WG 
have used rUp^ ~ 1.22 x lO^^GeV and pcri ~ S.l/i^lQ-^^GeV with h ~ 0.72 {h is the 
reduced present time Hubble parameter). If ruy:^ = leV, one obtains Qmin — 2.6 x 10~^^. 
Therefore, the minimum is located at tiny values of the quintessence vev (compared to the 
Planck mass). This makes a very important difference compared to the standard SUGRA 
case where, as already mentioned before, the minimum is close to mpj. Notice that, despite 
the fact that the vev of the quintessence field is now very small in comparison to the 
Planck scale, this does not mean that supergravity is no longer a necessary ingredient in 
the present context. This is because, when susy breaking is taking into account, the global 
susy limit is more subtle than simply taking the limit rUp^ +oo. 
Then, the scale M is given by 

4+a 

-] 



= (8vr)-"/^e,J ^ ) ^i8n)-'^/'Ql,^^. (4. 



PI 



If one uses the expression of Qmin deduced from Eq. ( [4.7| ) and the expression of the mass 
M, M"^"*"" ~ A^^^m", one obtains (assuming, as usual, A ~ 1) the expression of the scale 
^, namely 

.0 \ -2 



1/4 / X a/4 



^3/2 \ Peri 



a/8 

(4.9) 



This expression should be compared with Eq. (|l.3| ). We see that the difference lies in 
the presence of the term in the squared bracket which contains the gravitino mass and 
the critical energy density. This has disastrous consequences. If, as before, one chooses 
^3/2 ~ lOOGeV, then for, say, a = 4 one gets i/rrip^ ~ 10"'^*'. Even more serious, it is 
clear that there is no value of a which, for a reasonable value of the gravitino mass, would 
allow us to obtain a scale ^ greater than the TeV as it was the case for Eq. (|1.3|). Recall 
that the energy scale associated to the cosmological constant is A/nip^ ~ 10^^'^. This 
means that, in the case of quintessence, we have to build a complicated model and that, in 
addition, we have to fine tune the basic scale of the model even more than in the case of the 
cosmological constant where nothing else is needed. The origin of the problem can really 
be traced back to the fact that modeling quintessence in a more realistic fashion (breaking 
supersymmetry properly, taking into account the interaction between the various sectors 
etc ...) has modified the shape of the potential such that the original success of finding 



reasonable ^ from the "see-saw" formula (1.3) has vanished. 



Finally, the scale T can be expressed as 



^ ^ - , (4.10) 

PI 
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and this implies that "^^=1 ~ ^Z'^- This equation was used in Ref. [16]. 

The previous estimates rest on the assumption that, very quickly, the quintessence 
field is stabilized at its minimum. Considering to the drastic consequences evoked before, 
we now carefully check that this is indeed the case. 

4.2 Dynamics of the Quintessence Field 

Let us assume that the initial value of the field is such that Qua ^ Qmm- In this case, 
the potential studied above approximatively reduces to the Ratra-Peebles potential. In 
fact, it seems necessary to start from the branch Q~°' rather from the branch in order 
to have insensibility to the initial conditions. Indeed, if the potential is made of a series 
of monomial with positive powers, it is known that a fine-tuning of the initial conditions 
becomes necessary. In the case of the Ratra-Peebles potential, there is a particular attractor 
solution given by 

/ \ 3(1+^^3 )/("+2) 
Qattra = Qp iy— j , (4.11) 

where we have assumed that the field is a test field and that the background behaves as 

/ „ \ 2/(1+3^b) 

a(r?) = ap i^-^j , (4.12) 

r] being the conformal time and lo^ the equation of state parameter of the background fiuid 
(radiation or matter, that is to say lv^ = 1/3 or = 0). The time ryp is arbitrary and 
can be chosen at convenience. In practice, we will often consider that this is the time of 
reheating. The constant Qp can be expressed as 

18 l-u^l 

where denotes the Hubble parameter at time rj = r]p. loq is the equation of state 
parameter of the attractor and is given by ojq = (—2 + au}-^)/{a + 2). 

At this point, it is worth recalling the following point, already discussed at the end 
of the previous section. For = 0, the solution ( 4.12 ) a(r/) oc rf' is obtained from 



/Jmat oc \/a? . However, in the present context, one has to be more precise as in Eq. ( |3.25|) a 
dependence in the quintessence field is present. However, the Q dependence, at very small 
values of Q, is very weak. In other words nia (kQ^) cte if Q ^ m-p^. Therefore, in the 
regime under consideration, one can safely assume that pmat oc 1/a^ and the corresponding 
behavior of the scale factor follows. In other words, as already discussed, one can safely 
neglect the fact that the model is a chameleon in the present context. 

The attractor solution is completely specified once the fact that quintessence represents 
70% of the critical energy density today has been imposed. Let us evaluate its value just 
after infiation. At reheating, Zj-ch = 

10^^, the value of the field is 

Oattra (^.eh) ^ l^'^^' ' ^'^^'^ X Q'^J^^^^'^ , (4.14) 
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where, for simplicity, we have not considered the factor 9(1 —uJq)/ {20^) which is, roughly 
speaking, of order one (in the following we will always neglect this kind of factors since 
we are interested in order of magnitude estimates only). We have also taken the radiation 
contribution today to be ~ 10~^. For m^j^ = 100 GeV and a = 6, Qattra(-Zreh) ~ 
8.6 X 10~^^mpj. For rn^^i^ = leV and the same parameters, Qattra(-2rch) ~ 1-5 x lO^'^'^rripj. 
If we now compare the initial quintessence energy density with the energy density of the 
background (i.e. the energy density of radiation), one obtains 



f^g... M - 10-^°^/(^+^) X Q^^"^ . (4.15) 

For 7713^2 = 100 GeV and a = 6, one obtains ^Q^ttrS^reh) ~ 2.4 x 10~^^. For 771,3^2 = leV 
and a = 6, one obtains ^Q^tuS^teh) ~ 7.4 x 10"'''^. For comparison, we recall that, at 
reheating, one has Pg ~ 10^'^'^pcri- The previous estimates show that the energy density 
of the attractor is always initially very small in comparison with that of the background. 
This is not the case in the "standard" Ratra-Peebles case. This difference is due to the 
fact that the scale M now depends on Qmin which turns out to be a very small quantity. 



The solution given in Eq. ( 4.11 ) is valid as long as Qattra is small in comparison with 



Qmin and breaks down when Qattra ~ Qmm- This happens when 

2/(3+3a;„) 



^min / Hp 



(4.16) 



that is to say at a redshift given approximately by (strictly speaking, this expression is 
valid if Zjnin > 10^ since we used that lo^ = 1/3; otherwise, one has to take into account 
that the background becomes matter dominated) 

1 + .min ~ 10^^ X ~ 10 X O^f . (4.17) 

For our typical examples with m^j^ ~ 100 GeV, one obtains z^[^ — 2.4 x 10^3 to be 
compared with the reheating redshift, Zj-eh ~ 10^^. For rriy^ = leV, one gets Zmm — 7.7 x 
10^^. In both cases, the field reaches the minimum well before Big Bang Nucleosynthesis. 

Now, the field is of course not necessarily on the attractor initially. It is therefore 
important to estimate at which redshift the attractor is joined and to compare this redshift 
to 

^min- If Qini > Qattra (undcrshoot) then the field remains frozen. Therefore, the redshift 
at which the attractor is joined is given by the condition Qi^i = Qattra which results in 



Qundcr ( ^2p+2 



1/(3+3c^b) 



Qtirj^, ■ (4.18) 



Using this relation, one finds 

(- 7, ]„„ ~ 1 X I 

Q 



1 + Sunder ~ 10 X ( ^X'^ Q7^^ . (4.19) 



ini 
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One can check that, if the initial value of the field is the value on the attr actor given by 
Eq. (|4l^ ), then 

^under ~ ^reh as expected. More generally, we conclude that each time the 



field starts in a undershoot situation, that is to say initially Qattra < Qini < Qmin, it will 
join the attractor before reaching the minimum. 

Let us now consider the case of an overshoot, i.e. Qini < Qattra- Then, the field 
is first kinetic dominated until the potential energy becomes equal to the kinetic energy. 
When this happens, the field becomes frozen until the attractor is joined. During the phase 
dominated by the kinetic energy, the field behaves as 



Q = Qini + mp,^^(l-^) . (4.20) 

This allows us to estimate the redshift at which the field b and to 

compare it with Zmin- The redshift at which the field becomes frozen can be estimated to 
be 

1 + ~ 10^° X X n-Q^^^'^/' . (4.2f ) 

On the other hand, if the quintessence field behaves as in Eq. ( ^30|) , it reaches the minimum 
at a redshift given by 

1 + Zkin^min ~ 10^° X QqI^^ . (4.22) 

Therefore, the field becomes frozen before reaching the minimum if l + Zfroz > l + ^km— »mirn 

~ 1/2 

a condition which amounts to Qmin > ^q. .• As an illustration, let us consider again the 
model for which m-g^g — 100 GeV and a = 6. From the previous considerations, one knows 
that Qmin — 2.6 X 10~^^ and that, in order to have overshoot, 10~^^ < ^Qini ^ 10^^, this 
last bound corresponding to equipartition initially. Therefore, for the above condition to 
be satisfied, one needs iln ■ lO"^'^. 

' Vim ^ 

We conclude that, on more general grounds, the above condition cannot be satisfied 
unless the field is, at the beginning, very close to the attractor. Therefore, generically, 
when we have overshoot, the field reaches the minimum in the kinetic dominated phase 
and has no time to freeze out. 

4.3 Approaching the Minimum 

We now describe the behavior of the quintessence field when it starts feeling that the 
potential has developed a minimum. When the field is close to the minimum, the potential 
can be approximated by 



ViQ) ^ Klin + 2 [^3/2) iQ - QrmnY ■ (4.23) 

The mass of the field is the gravitino mass as established before and in Ref. [16]. Let us 
first consider the situation where the quintessence field is a test field. The Klein-Gordon 
equation, written with the number of e-folds as the time variable, reads 

d2 / Q \ / 1 dH\ d / Q \ ml^ df , , 



dN"^ \m^J \ HdNjdN\mpJ H^d{Q/mp^) 



-25- 



where V{Q) = rn^^f{Q). The Hubble parameter is given hj H = //p exp[— 3(1 + LL'g)A'^p/2], 
where Np is the number of e-folds counted from the time f] = rjp. The above equation can 
be easily integrated and the solution reads 



Q Qr 



m. 



^g3(l+a;B)iVp/2 



+A2J_(i_^ )/(2+2a;^) 



3(l + a;J Hp 



3/2 3(l+a;B)iVp/2 



(4.25) 



where Ju{z) is a Bcssel function of order u. As expected the field start oscillating around 
its minimum when its mass equals the Hubble parameter. Using the expression of the 
gravitino mass given before, one easily checks that this happens at a redshift of 

^lOx (g^in)"'/'. (4.26) 



When the mass is smaller than the Hubble parameter, the field is essentially frozen. Using 
the asymptotic expression of the Bessel functions for small arguments, one obtains 

Q Q m 

rsj 



+ Ai + A2e 



-3(1-0.3 )Afp/2 



(4.27) 



where Ai and A2 are two new constants, different from Ai and A2. Very rapidly, the 
branch proportional to A2 becomes negligible. In the opposite situation, i.e. when the 
mass is much larger than the Hubble parameter, one can use the asymptotic expansion of 
the Bessel functions for large values of the arguments and one obtains 



Q Q 



mm g-3(l- 



i)^p/2, 



Ai cos 



^ 3(1+0.3 )iVp/2 _ 7r(l -^b) _ ^ 
3(1 Hp 4(1+ a; J 4 



+A2 cos 



3(1 +0;^) Hp 



.!!^g3(l+^B)7Vp/2 ^ _ ^ 



4(1+0; J 4 



(4.28) 



where Ai and A2 are new constants. The oscillations are damped by a factor a~^^^~^B)/^ . 

Let us summarize the two possibilities. In case of an undershoot, the field joins the 
attractor and then reaches the minimum without any oscillatory phase If 
there is an overshoot, the field has no time to freeze out and goes directly from the kinetic 
dominated phase to the oscillatory phase. 

At some point, the quintessence field starts dominating the matter content of the 
Universe. In this case, the above treatment breaks down since the quintessence field is no 
longer a test field. Assuming, for simplicity, that p ~ Knin, the Klein-Gordon equation can 
still be solved explicitly. The solution reads 



Q_ 



+ e-3^/2<^Si 



cos 




m 



3/2 



+B2 sin 




m 



3/2 



StT (F„nn/m4)l/2 



N 



N 



(4.29) 
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and we still have very rapid oscillations, damped by a factor a~^/'^. For — 100 GeV, 
the dimensionless frequency {N being the time variable) is ~ 10'^^. Again the oscillations 
stop rapidly as their amplitude decays exponentially with the number of e-folds. 

4.4 Numerical Integration 

Of course, rather than the approximate considerations developed above, a full numerical 
integration would allow us to obtain the exact solution. Unfortunately, the realistic values 
of Qxain are so small that a simple Fortran code cannot handle the corresponding solution. 
However, one can check the previous analytical estimates for values of Qmin which are 
numerically reasonable (but not physically realistic). Then, having checked and validated 
the previous estimates, we will use them in a physically relevant situation. 

Let us consider a situation where the free parameters of the potential are given by 
Omin = 10"'^, a = 6 and B = \. This implies that the mass scale M ~ 8 x 10~^^mpj, the 
gravitino mass rn^i^ ~ 2 x 10~^^mpj and T ~ 2 x 10~'^^mp|. According to the previous 
estimates, the initial value of the field on the attractor is Qattra ~ 2.3 x lO^^^rripj corre- 
sponding to i^Qjittra ~ 1-8 ^ 10"'^^. On the attractor, the minimum of the potential is felt 
at a redshift of 1 + z^\^ ~ 1000. 

Let us now consider the initial conditions corresponding to equipartition, i.e. ^q-^^^ = 
10~^. This implies that Qini ~ 6.8 x 10~^^mpj and we have overshoot since Qini < Oattra- 
As a consequence, as explained before, the initial evolution is dominated by the kinetic 
energy and we have 

Q = Qini + "^p:/^§^ (1 - e-^) , (4.30) 

where is the total number of e-folds counted from reheating. Very quickly, we have 
Q = y^30Q;^;/(47r) ~ 0.00489mpj (or Q ~ 0.0245mpi). Let us notice that this value is 
greater than the value of the minimum. Although the field seems to be frozen, its time 
variation is still sufficient for the corresponding kinetic energy to be greater than the critical 
energy (a similar situation arises in the standard Ratra- Peebles scenario, see Ref. [14]). In 
fact, the field rolls down the potential so quickly that it goes through the minimum while 
the kinetic regime goes on (when the kinetic energy dominates, the fact that the potential 
has a minimum is irrelevant). The kinetic energy reaches the critical energy, by definition 
at Zkin^min) when the field is on the other side of the potential. As a consequence, the field 
will approach the minimum "from the right" . 

In the Ratra-Peebles case, the kinetic regime comes to an end at the redshift 1 + Zfroz ~ 
4.64 X 10^ while, in the present case where the potential possesses a minimum, this one 
is felt by the field at 

^kin— >min ~ 4.64 X lO^''. Therefore, in this case, the field does not 
enter the potential dominated regime at all and directly goes from the kinetic dominated 
regime to a regime where the minimum is felt and where the solution ( [1.25 ) is relevant. 



However, the solutions before -Zkin^min ~ 4.64 x 10^'^ given by Eq. ( [4.30| ) and after, given 
by Eq. ( 4.27| ) are the same. As a consequence, the coefficients Ai and A2 are such that 



the solution after Zkin^min is still given by Eq. ( [4.301) . So, even after Zkin^miiD the field 
remains "frozen" at Q ~ 0.00489mpj. Nevertheless, the evolution of the energy density 
changes and, instead of /> oc a~^, we have p ~ cte. 
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Figure 5: Evolution of the scalar field for the parameters, Qmin = 10^'', B = 1 and a = 6 with 
the initial conditions such that r^Q;,,; = 10~^ (equipartition and overshoot). The dashed horizontal 
represents the value Q = y^3r2Q~7/(47r). It is clear from this plot that the numerical estimates 
presented in the text are fully compatible with the numerical evolution. In particular, one can 
check that the oscillations start at z ~ 10^ and that, eventually, after rapid oscillations, the field 
stabilizes at its minimum. 



At 1 + 

^osci ~ lO'^, the Hubble parameter is equal to the gravitino mass and the damped 
oscillations of Eq. (|4.28 ) start. Then, the field stabilizes at its minimum. The evolution 
of the field is represented in Fig. |5| and the corresponding quintessence energy density is 
plotted in Fig. |. 

Let us now consider the case of an undershoot, i.e. the case where the initial value of 
the field is greater than the initial value on the attractor. We take rig.^. = 10~^^ which 
implies that Qini ~ 4.6 x lO^^^rripj. The corresponding numerical evolution is represented 
in Figs. 1^ and ^. Since we have undershoot, the field will be initially frozen until it reaches 
the attractor. This happens at Zunder ~ 4.6 x 10^^. On the attractor, the effective equation 
of state is ujq = (—2 + auj^)/{a + 2) = for a = 6 and = 1/3. As a consequence, 
the energy density of quintessence scales as matter, as can be checked on Fig. ^, until the 
minimum is reached. The minimum is felt at Zmin ~ 1000 which is also Zosci; the redshift 
at which the Hubble parameter is equal to the gravitino mass and the oscillations start. 
Therefore, we expect no phase where the field is frozen, as predicted by Eq. ( 4.27] ) but 
expect the oscillations to start immediately after the field has left the attractor. However, 
in the case of an undershoot, when the minimum is felt we necessarily have Q ~ Qmin- 
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Figure 6: Evolution of tlie quintessence energy density (solid line) for Qmin — 10^**, B — \ and 
a = 6 with the initial conditions such that J^Qj,,; = 10"'^. The dotted line represents the evolution of 
the radiation energy density while the dashed line is the matter energy density. The energy density 
freezes at Zkin^min and not at 2;froz as predicted in the text since the minimum of the potential 
is felt before the kinetic energy becomes equal to the potential energy. The analytical estimate 
Zkin^min ~ 5 X 10^" is in relatively good agreement with the actual value observed in the figure. 
The difference (about one order of magnitude) is probably due to the fact that the minimum is 
felt before the energy density is equal to pcri (which was the criterion used in order to derive the 
expression of ^kin^min) which has the effect to increase Zkin^min- 

In the previous case of an overshoot this was not the case because the kinetic energy of 
the field was dominating the potential energy at the moment where the presence of the 
minimum is seen by the field. Therefore, the amplitude of the oscillations is very small 
and, in practice, we see no oscillations at all. The field directly stabilizes at this minimum. 
This behavior is indeed observed in Fig. |^. 

As a conclusion of the subsection, we have checked that the numerical estimates derived 
before are confirmed by a numerical calculation in the (physically unrealistic) case where 
this one is possible. 

Let us summarize our findings in the realistic cases where the gravitino mass is of 
the eV or 100 GeV scales. In both cases the minimum of the potential is at extremely 
small values compared to the Planck scale. This implies that we can safely neglect the Q 
dependence of the particle masses and trust the supergravity expansion scheme in 1/mc. 
On the other hand, the existence of a minimum with a mass of the order of the gravitino 
mass leads to a drastic modification of the evolution of the quintessence field. In both 
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Figure 7: Evolution of tlie scalar field for the parameters, Qmin = 10~^, B = 1 and a = 6 with 
the initial conditions such that = f 0^"'^ (undershoot). The numerical estimates presented in 

the text are in good agreement with the numerical evolution. 

cases, either overshoot or undershoot, we find that the quintessence field settles down at 
the bottom of the potential before the beginning of BBN. Prom the point of view of late 
time physics, the quintessence sector plays the role of an effective cosmological constant 
with equation of state loq = — 1. Observationally this is not a problem since a cosmological 
constant is perfectly compatible with the currently available data. However, conceptually, 
we consider it as a second fundamental difficulty (the first one was he magnitude of ^) of 
the realistic model where supersymmetry breaking is taken into account. Indeed, it seems 
clear that the justification for building a complicated model which is just equivalent to a 
cosmological constant is very weak. One possible way out is to study whether the behavior 
of the perturbed quantities allows us to distinguish this model from a pure cosmological 
constant. In particular, although the equation of state is loq = —1, the mass of the 
quintessence field is now and, therefore, the corresponding Jeans length is very small 
compared to the usual case where it is the Hubble scale. This is why, in the next section, 
we analyze whether the perturbations of the quintessence field can give information on the 
dynamics of quintessence prior to BBN. 

4.5 Cosmological Perturbations 

In this section, we study how the perturbations of the quintessence field behave. In con- 
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Figure 8: Evolution of the quintessence energy density (solid line) for Qmin — 10^**, B — \ and 
a = 6 with the initial conditions such that ilgj^j = 10"'*^. The dotted line represents the evolution 
of the radiation energy density while the dashed line is the matter energy density. Initially, the 
energy density is frozen until the attractor is joined at Zundcr ^ 4.6 x 10^^. At 2:,nin = •Zosci, the 
presence of the minimum is felt by the field which stabilizes at its minimum. 



formal time, the perturbed Klein-Gordon equation Fourier space reads 



5Qk = 4QVk-2a2^$k, 



where 7i = a' /a and <I> is the Bardeen potential which describes the metric perturbations. 
We want to study the situation where the background field stands at its non-vanishing 
minimum, this field being a test field, i.e. the scale factor being given by the expres- 
sion (HH ). In this case, the right hand side of the above expression vanishes and the 
perturbed Klein-Gordon equation reduces to 



(l + 3a;Jr? 



^3/2 



6Qi, 



0. 



(4.32) 



It is clear that the Jeans mass of the field is now (^^3/2)"^' Since this scale is much smaller 
than the Hubble length, i.e. the Jeans length in the standard Ratra-Peebles case, this 
raises the question as whether quintessence could collapse and develop structures at very 
small scales. Unfortunately, the above equation cannot be solved explicitly. However, it 
can be analyzed in the case where the wavelength of the fluctuation is either larger or 
smaller than the Jeans length. 
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Let us start by assuming that fcph ^ ^-3/2' wavelength of the Fourier mode is 



smaller than the Jeans length. Then the term proportional to (771,3^2)^'^^ perturbed 
Klein-Gordon equation can be neglected. In this case, the solution reads 



5Qk = Ai{k)r]''J^{kr]) + A2{k)rj''J^^{krj) 



(4.33) 



where Ai{k) and A2{k) are two integration constants fixed by the initial conditions which 
are not important for us in the present context and where 3^ is a Bessel function of order 
u which is a function of the background equation of state only 



3fa;^ 



11 



2(1 + 3^;^ 



(4.34) 



Moreover, kph ^ ^3/2 implies fcph ^ H and then using the asymptotic behavior of Bessel 
functions for large values of their argument, one obtains the approximate expression 



Ai{k) cos [krj - ^ - + A2{k) cos (jirj + ~ ^ 



(4.35) 



We have oscillations since we consider modes of wavelength of smaller than the Jeans 
length. The overall amplitude behaves as rf~'^l'^. Since one has 



l + Sw, 



<0, 



(4.36) 



there is no growing mode in this situation. This result makes sense as the pressure forces 
counter balance the gravitational force which tends to make the system collapse. 

Let us now consider the situation where /cph ^ "^3/2' ^-^^ where the wavelength of the 



3/2' 

Fourier mode under consideration is larger than the Jeans length 



m 



3/2^ 



In this case. 



ignoring the term A;^ in Eq. ( [4.32 ), the equation for bQ\^ can also be integrated exactly. 
It is interesting to compare this equation with the standard equation (i.e. in the case 
where the potential is the Ratra-Peebles one, V = M^"'""Q~") for JQk on large scales 
when the field is on the attractor. This case was studied in Ref. [14]. On the attr actor, 
the second derivative of the potential is given by d^F/dQ^ = 9-?/^(a + 1)(1 — lJq)I(2q)^ 
wq being a constant, and, as a result, the terms a^d'^V/dQ^ scales as r/~^. This is why the 
homogeneous equation admits simple power-law solutions. Here, the term d^V/dQ^ is a 
constant equal to rriy^ and, therefore the time dependence of the term a'^d'^V/dQ'^ is now 
given by a^, i.e. a power-law of the conformal time, namely r/^ for the radiation dominated 
era and r/^ for the matter dominated epoch. Therefore, we expect the solutions for dQ]^ on 
large scales to be different from the Ratra-Peebles case. Indeed, the result reads 



SQk = Bi{k)r] '+'-B J^ 
+B2{k)r] J_ 



l + 3c^B 
3(1 +Wb. 
1 + 3a; 



m 



3/2^P 



l + 3a; 



3{l+iJp.) 



B l+SiJj 



3(1+0; ) 



3(1 



(4.37) 
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where Bi{k) and B2{k) are two arbitrary constants and where the order ji of the Bessel 
functions can be written as 

(4.38) 



The same expression can be also rewritten in a form for which the physical interpretation 
is easier 



3(u>B-l) 



m 



■ 

3/2 



3(1 +0^3) H 



3(a;„-l) 



+ B2{k)r] 1+^"B J_ 



■ 
3/2 



3(1 +0^3) H 



(4.39) 



where rn^i^ S> H. Therefore, in the above equation, one can use the asymptotic expression 
of the Bcsscl functions for large values of their argument. This means that the overall 
amplitude behave as 77^K-i)/(i+3'^b)-V2_ xhis gives t?-^/^ ^.^^ ^-3/2 f^^, ^ 0,1/3. 

Again this function is a decreasing function of r]. Hence no growing mode is generated by 
perturbation of the quintessence scalar field. 

Together with the fact that for most of the history of the Universe, the quintessence field 
sits at the minimum of its potential, this implies that the dynamics of quintessence when 
coupling to supersymmetry breaking in a hidden sector leads to an absence of deviations 
from a pure cosmological constant in the late Universe, even at the perturbed level and 
even though the Jeans length of the field is now considerably smaller than in the standard 
Ratra-Peebles case. 



5. Conclusions 



We have presented a cosmological analysis of quintessence models in supergravity coupled 
to matter. For the quintessence sector, we have used the SUGRA model whose main 
feature is to reduce to the Ratra-Peebles potential at small values of the quintessence field. 
This is effectively the only property which is required here. Hence our results generalize 
to quintessence models with the same type of potentials. Supersymmetry is broken in a 
hidden sector which couples gravitationally both to the observable sector, i.e. the MSSM, 
and the quintessence sector. Requiring that the gravitino mass is large enough to lead to 
acceptable masses for the sparticles implies that the quintessence potential is drastically 
modified by the presence of the hidden sector. In particular, as soon as the hidden sector 
fields are stabilized and supersymmetry is broken, we find that the quintessence potential 
develops a minimum acting as an attractor for the quintessence field in the very early 
universe. Indeed the quintessence field settles down at the minimum of the potential before 
BBN and both at the background and perturbation level, the model becomes equivalent 
to a pure cosmological constant scenario. This is the first quintessential difficulty as no 
observational consequences seems to spring from such a scenario. Moreover, the energy 
scales of the quintessence sector required to obtain the correct vacuum energy now are 
minute as a consequence of the tiny value of the quintessence field now. As a result, one 
must introduce a highly fine-tuned scale in the quintessence sector which is as difficult to 
explain as the smallness of the vacuum energy. On the other hand, one may hope to preserve 
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the SUGRA potential runaway shape by tuning the hidden sector dynamics. In this case, 
we find that the local test of gravity (fifth force, weak equivalence principle and proton to 
electron mass ratio variation) are incompatible with the vev of the quintessence field implied 
by the presence of an attractor which, again, is required if we want insensitivity to the initial 
conditions. Hence it seems difficult to build models of quintessence in supergravity where 
the cosmological, gravitational and particle physics aspects are compatible. One of the 
plausible possibilities consists in getting rid of the regularity of the Kahler potential for 
small values of the quintessence field. This is what happens in no-scale models for instance. 
The analysis of no-scale models leading to both supersymmetry breaking and quintessence 
is left for future work [17]. 
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